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A POPULATION OF RELIC INTERMEDIATE-MASS BLACK HOLES IN THE HALO OF THE MILKY WAY 
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ABSTRACT 

If "seed" central black holes were common in the subgalactic building blocks that merged to form 
present-day massive galaxies, then relic intermediate-mass black holes (IMBHs) should be present 
in the Galactic bulge and halo. We use a particle tagging technique to dynamically populate the 
N-body Via Lactea II high- resolution simulation with black holes, and assess the size, properties, and 
detectability of the leftover population. The method assigns a black hole to the most tightly bound 
central particle of each subhalo at infall according to an extrapolation of the Mbh — 0* relation, and 
self-consistently follows the accretion and disruption of Milky Way progenitor dwarfs and their holes 
in a cosmological "live" host from high redshift to today. We show that, depending on the minimum 
stellar velocity dispersion, er m , below which central black holes are assumed to be increasingly rare, 
as many as ~ 2000 (a m — 3 kms -1 ) or as few as ~ 70 (a m = 12 kms -1 ) IMBHs may be left 
wandering in the halo of the Milky Way today. The fraction of IMBHs kicked out of their host 
by gravitational recoil is < 20%. We identify two main Galactic subpopulations, "naked" IMBHs, 
whose host subhalos were totally destroyed after infall, and "clothed" IMBHs residing in dark matter 
satellites that survived tidal stripping. Naked IMBHs typically constitute 40-50% of the total and are 
more centrally concentrated. We show that, in the a m = 12 kms -1 scenario, the clusters of tightly 
bound stars that should accompany naked IMBHs would be fainter than my = 16 mag, spatially 
resolvable, and have proper motions of 0.1-10 milliarcsec per year. Their detection may provide an 
observational tool to constrain the formation history of massive black holes in the early Universe. 
Subject headings: black hole physics - Galaxy: halo - galaxies: evolution - dwarfs - method: numerical 
- stellar dynamics 



1. INTRODUCTION 

Direct dynamical measurements show that most lo- 
cal massive galaxies host a quiescent massive black hole 
(MBH) in their nuclei. Their masses have been found to 
correlate tightly with the mass (Haring & Rix 2004) and 
the stellar velocity dispersion of the host stellar bulge, as 
manifested in the Mbh — a * relation of spheroids (Fer- 
rarese & Merritt 2000; Gebhardt et al. 2000). It is not 
yet understood whether such scaling relations were set in 
primordial structures and maintained throughout cosmic 
time with a small dispersion, or indeed which physical 
processes established such correlations in the first place. 
It is also unclear whether there exists a minimum host 
galaxy mass or velocity dispersion below which MBHs 
are unable to form or grow. The subset of the MBHs 
that populates the centers of dwarf galaxies and very late 
spirals and is undergoing accretion can be detected as ac- 
tive galactic nuclei (AGNs). Observations of AGNs with 
estimated black hole masses in the range 10 5 -10 6 M ap- 
pear to be consistent with an extrapolation of the local 
-Mbh — <?* relation for inactive galaxies down to stellar 
velocity dispersions of 30 kms -1 (Barth et al. 2005; Xiao 
et al. 2011). The existence of "intermediate-mass" black 
holes (IMBHs) with 20 M < M BH < 10 5 M Q remains 
in dispute. Dynamical mass measurements have shown 
the presence of IMBH candidates in the cores of the glob- 
ular clusters Gl (M BH = 1.8±0.5 x 10 4 M Q , Gebhardt et 
al. 2005), w Centauri (M BH = 4.7±1.0 x 10 4 M Q , Noyola 
et al. 2010), NGC 1904 and 6266 (M BH = 3±1 x 10 3 M 
and M BH = 2 ± 1 x 10 3 M©, Liitzgendorf et al. 2012). 
The X-ray spectra and bolometric luminosities of the ul- 



traluminous off-nuclear X-ray sources detected in nearby 
galaxies may imply the presence of IMBHs with Mbh <; 
500 M Q (Farrell et al. 2009; Kaaret et al. 2001). 

While the "seeds" of the MBHs powering the z > 6 
Sloan Digital Sky Survey (SDSS) very luminous, rare 
quasars (Fan 2006) must have appeared at very high 
redshifts and grown rapidly to more than 10 9 M Q in 
less than a Gyr, the problem of their origin and oc- 
cupation fraction in early galaxies remains unsolved. 
MBHs may have grown from the remnants of Popula- 
tion III (Pop III) star formation in sub-dwarf galaxies at 
z > 15 (e.g. Madau & Rees 2001; Volonteri et al. 2003; 
Tanaka & Haiman 2009), or from more massive precur- 
sors formed by the "direct collapse" of large amounts of 
gas in dwarf galaxy systems at later times (Loeb & Ra- 
sio 1994; Koushiappas et al. 2004; Begelman et al. 2006; 
Lodato & Natarajan 2006; Mayer et al. 2010). Mas- 
sive seeds have larger masses than Pop III remnants, but 
form in rarer hosts. Questions remain in the Pop III 
remnant model about the ability of ~ 100 M Q seed holes 
to grow at the Eddington rate for tens of e-folding times 
unimpeded by feedback, and in the direct collapse model 
about the needed large supply of low angular momentum 
gas that must accumulate in the center of young galax- 
ies before fragmentation and star formation sets in (see, 
e.g., Haiman 2012 and references therein). 

In this Paper we assess the size, properties, and de- 
tectability of the leftover population of IMBHtQ that 

1 In the following, we will use the term IMBH to refer to any 
hole with 20 Mq < Mbh < 10 5 Mq, i.e. more massive than the 
stellar-mass black holes found in X-ray emitting binary systems 
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is predicted to survive today in the halo of the Milky 
Way (MW) galaxy by two fiducial seeding scenarios. It 
was first pointed out by Madau & Rees (2001) that, if 
seed holes were indeed common in the subgalactic build- 
ing blocks that merged to form the present-day massive 
galaxies, then a numerous population of relic IMBHs 
should be present in the Galactic bulge and halo (sec 
also Volonteri et al. 2003; Islam et al. 2003; Volonteri 
& Perna 2005). Micic et al. (2006, 2011) used collision- 
less simulations to study the effect of gravitational re- 
coil kicks on the IMBH distribution in present-day galax- 
ies and test merger-driven recipes for black hole growth, 
van Wassenhove et al. (2010) ran Monte Carlo realiza- 
tions of the merger history of massive galaxy halos to 
study IMBHs in MW sat ellites at z = 0. A sim ilar ap- 
proach was employed by lO'Learv fc Loebl (pOOl to es- 
timate the expected number of recoiled IMBH remnants 
present today in the MW halo. Cosmological hydrody- 
namic simulations of the growth of MBHs are either lim- 
ited by resolution to galaxy hosts with Mhaio > 10 10 M Q 
(e.g. Di Matteo et al. 2008) or stop at high redshift (e.g. 
Bellovary et al. 2011). In all approaches seed holes are 
planted by hand. To improve on and complement the 
above calculations, we use here a particle tagging tech- 
nique to dynamically populate the N-body Via Lactea II 
(VLII) extreme-resolution simulation with IMBHs. As 
we shall discuss, this method allows us to self-consistently 
follow the accretion and disruption of MW progenitor 
dwarfs and the kinematics of their holes in a cosmologi- 
cal "live" host, and the build-up of a large population of 
Galactic "naked" wandering IMBHs. 

2. BLACK HOLE TAGGING TECHNIQUE 

The cosmological ACDM VLII simulation, one of the 
highest-precision N-body calculations of the assembly of 
the Galactic halo to date (Diemand et al. 2008), was 
performed with the PKDGRAV tree-code (Stadel 2001). 
It employs just over one billion 4, 100 M particles to 
model the formation of a M 2 oo = 1.93 X 10 12 M Q Milky 
Way-sized halo and its substructure. About 20,000 sur- 
viving subhalos of masses above 10 6 M Q are resolved 
today within the main host's r2oo = 402 kpc (the ra- 
dius enclosing an average density 200 times the mean 
matter value). Central black holes are added to subha- 
los following the par ticle tagging technique detailed in 
IRashkov et al.1 (|2012t ) and quickly summarized here. In 
each of 27 (out of the 400 available) snapshots of the sim- 
ulation, chosen to span the assembly history of the host 
between redshift z = 27.54 and the present, all subhalos 
are identified and linked from snapshot to snapshot to 
their most massive progenitor: the subhalo tracks built 
in this way contain all the time-dependent structural in- 
formation necessary for our study. We then: 1) identify 
the simulation snapshot in which each subhalo reaches 
its maximum mass, Mhaio, before being accreted by the 
main host and tidally stripped; 2) measure the subhalo 
maximum circular velocity V max ; 3) link it to the stel- 
lar line-of-sight velocity disp ersion, cr*, using the rela- 
tion Vmax = 2.2(7, derived by IRashkov et al.l (|2012fl ; and 
4) tag the most tightly bound central particle as a black 
hole of mass Mbh according to an extrapolation of the 

(Orosz et al. 2007) and 40 times less massive than Sgr A*. 



Mbh — relation of Tremaine et al. (2002), 

^ = /l0 6 ' 91 (loo^) 4 (-^kms- 1 ) 
M Q 1 100 (er* < 6 kms- 1 ). v ' 

Note that significantly steeper Mbh — c» relations have 
been derived recently by, e.g. Graham et al. (2011) and 
McConnell & Ma (2013). We will discuss the impact of 
a steeper power-law on our results in § 6. 

By neglecting all subhalos with Mhaio < 10 7 M Q , we 
restrict our analysis to 3,204 such tracks. Any evolution 
of the tagged holes after infall is purely kinematical in 
character, as their satellite hosts are accreted and dis- 
rupted in an evolving Milky Way-sized halo. After tag- 
ging, the black hole particles are tracked down to the 
z = snapshot. The main host is assigned a central 
black hole at z — of mass equal to that of Sgr A*, 
M BH = 4 x 10 6 M Q (Chez et al. 2008). While we tag at 
most one hole per subhalo, we allow multiple systems to 
form when subhalos merge after infall (see below) . 

Figure Q] shows the distributions of stellar velocity dis- 
persion at infall (left panel) and time since infall (right 
panel) for the VLII subhalo population. The color cod- 
ing separates self-gravitating subhalos that survive their 
accretion event and become MW satellites from those 
that are totally disrupted. The latter are found to fall in 
preferentially at earlier times, and have a median infall 
redshift of 4.5. 

3. DEMOGRAPHY OF IMBHS 

Below we discuss two simple models of seed hole for- 
mation that may be illustrative of more realistic growth 
scenarios, and differ only in the black hole occupation 
fraction as a function of the host's stellar velocity dis- 
persion. 

3.1. Population III Remnants 

Let us first consider a scenario where the black hole 
occupation fraction is of order unity at infall in all sub- 
halos with stellar velocity dispersions > a m — 3 kms -1 
and drops to zero below a m . This value is compara- 
ble to the stellar velocity dispersion measured today in 
the ultra-faint MW satellite Segue 1 (Simon et al. 2011), 
i.e. this model places seed holes in small-mass subhalos 
that are know n to have been rathe r inefficient at form- 
ing stars (e.g. R ashkov et al.ll2012t Cooper et al. 2010; 
Koposov et al. 2009). The mass distribution of subhalos 
with <7* > 3 kms -1 at infall has a median of 5 x 10 7 M Q . 
We assign black holes masses following the relation ([1]). 
As most of the subhalos in Fig. [T] host then a 100 M Q 
black hole, we refer to this model as "Pop III remnants" . 

Two large subpopulations of Galactic black holes can 
then be readily identified today within r2oo , 798 "naked" 
IMBHs, whose host subhalos have been totally disrupted 
after infall, and 1,234 "clothed" IMBHs residing in satel- 
lites that survived tidal stripping. Within the latter cat- 
egory, 1,096 are single black hole systems, i.e. are located 
in subhalos that host only one hole at the present epoch. 
The others arc multiple systems: we find 42 dwarf subha- 
los with 2 holes, 7 with 3, 5 with 4, and 1 with each of 6 
and 7 holes. These multiple systems are produced when 
their host subhalos merge together after infall, as they 
fall into the main halo in groups. We assume that by 



Rashkov & Madau 



3 




10 20 30 40 4 8 12 

cr* [km/s] time since infall [Gyr] 

Fig. 1. — Histograms of the stellar velocity dispersion at infall (left panel) and infall lookback time (right panel) for all subhalos with 
A^halo > 10 7 Mq in the VLII simulation. The color coding separates self-gravitating subhalos that survive their accretion event and become 
MW satellites (blue) from those that are totally disrupted (red). The vertical lines in the left panel mark the minimum velocity dispersion 
below which central black holes are assumed to be increasingly rare in the two seeding scenarios discussed in the text. 



that time, ram pressure stripping has removed any left- 
over gas from the system, making it difficult for IMBHs 
to sink to the center of the potential well and merge. 
Substructures that host 5 or more IMBHs are found to 
be more than 250 kpc from the Galactic center and have 
masses in the range 2 x 10 9 — 3 x I0 9 M Q today. Figure 
[2] (left panels) shows an image of the projected spatial 
distribution of IMBHs in today's inner and outer Galac- 
tic halo according to this model, while Figure [3] depicts 
their mass function and radial distribution. Naked holes 
are more concentrated towards the inner halo regions as 
a consequence of the tidal disruption of infalling satel- 
lites. Indeed, within 10 kpc of the center, most IMBHs 
are naked. The second most massive hole (after "Sgr 
A*") has M B h = 8.6 x 10 4 M (i.e. was tagged at in- 
fall in a large satellite with cr* — 32 kms -1 ), is naked, 
and is located at a distance of 55 kpc from the center. 
There are 48 IMBHs above M BH = 10 3 M Q , of which 
29 are naked. Most holes are assigned the seed mass of 
100 M©, i.e. they are not allowed to grow in this model: 
we count a total of 1,801 such "light" IMBHs, of which 
695 are naked. 

3.2. Direct Collapse Precursors 

Consider by contrast a scenario where more massive 
seeds form in rarer hosts. Models in which seed black 
holes become very rare for halos with -Mi la i < 10 9 M© 
have been recently discussed by, e.g., Bellovary et al. 
(2011) and Devecchi et al. (2012). Here, we take the hole 
occupation fraction to be of order unity in all subhalos 
with stellar velocity dispersions > a m = 12 kms -1 , and 
to drop to zero below a m . Black hole masses are again 
assigned following equation ((T|) down to cr m . The mass 
distribution of subahalos hosting central black holes at 
infall now has a median of 3 x 10 9 M Q , and the minimum 
hole mass is M B h ~ 1, 700 M©. 

This tagging prescription gives origin to 39 naked and 
33 clothed IMBHs (of which 27 are single systems). Fig- 
ure [5] (right panels) shows an image of the projected spa- 
tial distribution of IMBHs in today's inner and outer 
Galactic halo according to this model, while Figure [3] 



(left panel) depicts their mass function. There are 20 
IMBHs above M BH = 10 4 M Q , of which 13 are naked. 

4. BLACK HOLE RECOILS 

The asymmetric emission of gravitational waves pro- 
duced during the coalescence of a binary black hole sys- 
tem imparts a velocity kick to the system that can dis- 
place the hole from the center of its host. The magnitude 
of the recoil will depend on the binary ma ss ratio and 
the d irection and magnitude of their spins (| Baker et al.l 
120081 ) , but not on the total mass of the binary. When the 
kick velocity is larger than the escape speed of the host 
halo, a hole may be ejected into intergalactic space before 
becoming a Galactic IMBH. The demography of IMBHs 
discussed in the previous section accounts for recoiling 
holes as follows. 

Following the results of high-resolution self-consistent 
gasdynamic al simulations of binary mergers of disk 
galaxies by iKazantzid is et al.l (|2005h . we assume that 
unequal-mass galaxy mergers with mass ratios larger 
than 4:1 will not lead to the formation of a close black 
hole pair at the center of the remnant, and hence of a 
recoiling hole. All major mergers (mass ratios smaller 
than 4:1) produce instead recoiling holes with a kick 
velocity that is equal to the most likely value of the 
probability distrib ution function for r andomly oriented 
spins (see Fig. 3 of lGuedes et al.ll201lD . In practice, this 
means that all recoiling holes will have speeds in excess of 
100 km s _1 and will be ejected from the system. This is a 
conservative assumption in the sense that it tends to min- 
imize the actual number of relic IMBHs in the Galactic 
halo. We trace all <4:1 halo mergers following the VLII 
merger tree back to redshift 10. This was built by identi- 
fying a ll subhalos in each sn apshot with the 6DFOF halo 
finder (jDiemand et al.1 12006). and linking them from one 
snapshot to the next by the id numbers of the shared 
particles. We identify all halos that would have been as- 
signed a hole at infall, but had in their past a <4:1 major 
merger with another black hole-hosting system. As the 
kick velocity is always larger than the escape speed from 
the host, the hole is removed from the final catalog of 
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Fig. 2. — Projected distribution of IMBHs in the Galactic halo today. Left panels: Pop III remnant model. Right panels: Direct collapse 
progenitor model. The color coding is as in Fig. [T] while the size of the points is proportional to log 10 Mbu- The stellar halo from Rashkov 
et al. (2012) is plotted in the background. Top row: 200 kpc box. Bottom row: Zoom-in of the inner 20 kpc region. 



Galactic IMBHs. 

The size of the recoiling population is clearly very sen- 
sitive to the hole occupation fraction. Figure @] shows the 
total number of halo major mergers that produce recoil- 
ing IMBHs in the "Pop III" and "direct collapse" seed 
scenarios. The former produces 577 recoiling holes (22% 
of the total black hole population), mostly belonging to 
~ 10 8 M Q subhalos at infall. The latter produces only 5 
recoiling holes (6% of the total) all from rather massive 
hosts. Note that it is because of recoils that the num- 
ber of Galactic IMBHs heavier than a few thousand solar 
masses is actually larger in the direct collapse model (see 
Fig. [3]). The demography of Galactic IMBHs in our two 
scenarios is summarized in Tabic [1] 



5. SPATIAL DISTRIBUTION OF IMBHS 

The spatial distribution of Galactic IMBHs today is 
shown in Figure [5j Naked holes dominate in the inner ~ 
50 kpc from the halo center. Their number density, how- 
ever, plummets at larger distances, where most IMBHs 
are still hosted in surviving substructures. This is a man- 
ifestation of the fact that subhalos orbiting closer to the 
dense central regions of the main host are more likely 
to be tidally disrupted, and leave their IMBHs exposed. 
In the outer halo regions, the lower background density 
makes substructures more resilient to tidal forces, allow- 
ing them to survive and retain their holes. A comparison 
of the IMBH number density to the density profile of the 
dark matter shows that the cumulative distribution of 
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Fig. 3. — Left panel: Mass function of relic IMBHs within r2oo m the two scenarios discussed in the text. Solid histogram: Pop III 
remnants. Hatched histogram: Direct collapse precursors. Right panel: Location of "naked" (red diamonds) and "clothed" (blue circles) 
IMBHs in the AfpsH-Galactocentric distance plane for the Pop III model. Naked IMBHs are more concentrated towards the inner halo 
regions, a trend that is also observed in the direct collapse, massive seed scenario. 



TABLE 1 
Demography of Galactic IMBHs 

Model a m mcd(M halo ) Awaked ^clothed -^massive 

min(A^BH) max(AiBH) -^recoiled 

[kms-'l [M Q ] (A/ BH > 5,000 M ) [M ] [Mq] 

Population III 3 5 x 10 7 798 1234 21 100 1.9 x 10 5 577 
Direct Collapse 12 3 x 10 9 39 33 37 1691 1.9 x 10 5 5_ 

Note. — Column 1 indicates the seeding scenario, columns 2, 3, 4, 5, 6, 7, 8, and 9 give the minimum subhalo line-of-sight stellar velocity 
dispersion required for hosting a central black hole, the median halo mass at infall, the number of "naked" and "clothed" Galactic IMBHs, 
the number of Galactic IMBHs more massive than 5000 Mq, the mass of the lightest and heaviest Galactic IMBHs, and the number of 
recoiled IMBHs, respectively. 



Galactic IMBHs (naked plus clothed) follows the dark 
matter profile rather well, but that each individual sub- 
population does not. The inset in Figure [5] shows that 
the IMBH mass density profile does not follow the den- 
sity profile of the dark matter, as most of the mass in 
black holes is concentrated in the inner regions. 

Since naked IMBHs more massive than dark matter 
particles do not experience proper dynamical friction in 
our simulation, we have estimated a posteriori the effect 
of black hole orbital decay in two different ways - using 
the slowly-decaying circular orbit approximation, as well 
as computing t he instantaneous decelera tion from dy- 
namical friction (|Binnev &: T rcmaine 2008) in each avail- 
able VLII snapshot. We find that, even for the heavi- 
est IMBHs in the highest density regions of our simu- 
lation, dynamical friction timescales exceed the Hubble 
time at all poin ts along the orbit . A si milar conclusion 
was reached by lO'Learv fc Loeb ( 2009) in their study. 
Of course, our N-body simulation does not account for 
possible hydrodynamical effects and interactions with the 
galaxy stellar disk in the very center of the VLII halo. 

6. DISCUSSION 

We have used a particle tagging technique to populate 
the subgalactic building blocks of a present-day MW- 
sized galaxy with black holes, and assess the size, proper- 
ties, and detectability of the leftover population. Our ap- 
proach combines a computationally expensive, extreme- 
resolution, N-body simulation like VLII, in which struc- 
tures grow ab initio in a fully cosmological framework, 



with simple prescriptions for seeding, in post-processing, 
progenitor subhalos with IMBHs. Insofar baryonic mate- 
rial does not appreciably perturb the collisionless dynam- 
ics of the N-body component, the dynamical association 
of black holes with individual particles in the N-body 
component should correctly reproduce the spatial and 
kinematic properties of IMBHs in galaxy halos, in par- 
ticular of the naked subpopulation whose host satellite 
galaxies have been totally destroyed after infall. This 
level of detail is unavailable to a standard merger-tree 
approach. 

We have discussed two simple models of seed hole for- 
mation that may be illustrative of more realistic growth 
scenarios, a "Pop III remnant" model in which small- 
mass black holes populate subhalos with stellar velocity 
dispersion as low as 3 kms" 1 , and a "direct collapse pre- 
cursor" model in which holes become very rare in systems 
with stellar velocity dispersion below 12 kms -1 , and are 
more massive. It is important to keep in mind that the 
Pop III route does not necessarily require the formation 
of extremely massive stars from gas of primordial com- 
position. Already at metallicities Z < 0.01 Zq, the mass 
loss rates of massive stars through radiatively driven stel- 
lar winds are predicted to be rather small (with rates 
decreasing with metallicity asMoc Z 69 , see Vink et al. 
2001). If low-metallicity stars above 50 M Q collapse to 
black holes after losing only a small fraction of their ini- 
tial mass, then IMBHs with masses above the 5 — 20 M 
range of known "stellar-mass" holes may be the in- 
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Fig. 4. — Top panel: Number of subhalo major mergers that 
produce a recoiling IMBH, as a function of redshift. Light green 
curve: Pop III remnants model. Dark green curve: Direct collapse 
precursors. Lower panel: Mass distribution (at infall) of all the 
subhalos whose central black holes were ejected following a major 
merger. The color scheme is the same as in the top panel. 

evitable end product of early star formation. A standard 
Salpeter initial mass function (IMF), dN/dm* oc m* 2 
in the range 1 — 100 M Q , would produce an IMBH ev- 
ery 1000 M Q of stars formed. A top-heavy IMF like that 
found in recent simulations of primordial star formation 
by Stacy & Bromm (2012), dN/dm r oc m^ 017 , over the 
same 1 — 100 M© mass range, would generate an IMBH 
every 100 M Q of stars formed. Note that 1000 M of 
stars in a 5 x 10 7 M Q subhalo imply a star formation 
efficiency that is much too high to match the observed 
dearth of faint Milky Way satellites (e.g. Rashkov et 
al. 2012), and a top-heavy IMF must then be invoked for 
consistency. 

Galactic IMBHs may light up if they pass through 
dense regions of the galaxy and accrete from the in- 
terstellar medium (Volonteri & Perna 2005). IMBHs 
present in the lensing galaxy of multiply-imaged back- 
ground QSOs will produce monopole- or dipole-like dis- 
tortions in the surface brightness of the QSO images and 
may be detectable by next generation sub-millimeter tele- 
scopes with high angular resolution (Inoue et al. 2013). 
It is interesting at this stage to discuss the observability 
of the stellar cusps that may accompany naked IMBHs 
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Fig. 5. — Spatial distribution of Galactic IMBHs versus the dark 
matter density profile in VLII. Main panel: IMBH number density 
Inset: IMBH mass density. The mass density of the dark matter 
has been multiplied by a factor of 1,000 for comparison. Light 
green: Pop III remnants scenario, naked (solid curves) and clothed 
(dotted curves) IMBHs. Dark green: direct collapse progenitors 
scenario, naked (solid curves) and clothed (dotted curves) IMBHs. 
Gray: VLII dark matter density profile. The total number density 
of iMBHs follows the dark matter density profile, while their mass 
density does not. 
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Fig. 6. — Distribution of V-band peak surface brightness (within 
the central 1 arcsec) of stellar cusps around naked IMBHs in the 
"direct collapse" model. Only 20% of the initial stellar mass (equal 
to 2 X Mbh) is assumed to be retained today. To account for the 
partial SDSS sky coverage, the expected numbers should be scaled 
down by a factor of ^5. 



wandering in the Milky Way halo. After the complete 
disruption of their host subhalo, naked IMBHs are still 
expected to carry w ith them a Bahcall-Wolf stellar cusp 
(jMerritt et al.ll2009f ) which, upon formation and in the 
absence of a kick from gravitational recoil, has a mass 
of order the mass of the hole and an extent of order the 
hole's radius of influence, 



GMbh/o* 



(2) 



Throughout their dynamical evolution, these stellar clus- 
ters will lose mass by relaxation and expansion, as well as 
tidal disruptions by the black hole (|Komossa fc Merritt] 



Rashkov & Madau 



7 



120081 10'Learv fc Loebll2012l ). The relaxation timescale, 



tr w 10 9 



M. 



10 5 M, 



5/4 



yr. 



(3) 



ranges from several Myr for the smallest IMBHs to a Gyr 
for the largest. As IMBHs get accreted by the Milky Way 
early in the history of the Universe, the above processes 
can cause a 60-80% fractional loss o f stellar mass by the 
present day (|Q'Learv fc Loebll2012ft . At the same time, 
the physical extent of the clusters is expected to grow like 
t 1 / 3 a fter the first relaxation timescale ()0'Learv k, Loebl 
12012ft . resulting in a total expansion by a factor of 16 (for 
M. ~ 10 3 M Q ) to 3 (for M. - 5 x 10 4 M Q ) over a Hub- 
ble time. Assuming that the initial extent of the clusters 
is dictated by the IMBH's sphere of influence at infall, 
their physical extent today would range between 0.5pc 
and l.Opc. Despite their expansion, stellar clusters sur- 
rounding IMBHs would still be very compact today and 
may appear point-like or extended depending on their 
distance. Photometrically, they would have colors similar 
to those of old stars, while spectroscopically they would 
be distinguishable by the high velocity dispersions. 

We have modelled the apparent brightness distribu- 
tion of stellar cusps around naked IMBHs in the di- 
rect collapse scenario, which produces the largest num- 
ber of massive IMBHs today. We assume an initial stel- 
lar mass that is twice the mass of the black hole and a 
present-day stellar mass-to-light ratio of 4 (as appropri- 
ate for_a^i_old_s^lai_population with a Salpeter IMF - 
see lRashkov et afll2012L for details). The mass density 
in the star clusters is expected to follow a steep power 
law with drop-off power a ~ -2.15. Given these physical 
characteristics, we calculate the peak surface brightness 
(within a typical SDSS seeing element of angular size 
~ 1 arcsec) of the clusters as would be seen by a mock 
observer located 8 kpc from the center of the VLII main 
halo. The distribution is shown in Figure[B] The angular 
size out to which the brightness profile will be brighter 
than 25 mag arcsec" 2 is shown in Figure [3 Most clusters 
would be resolvable, as their angular extent spans 2 to 
10 1-arcsec resolution elements. In these calculations, we 
have accounted for the fact that the SDSS points away 
from the Galactic Center by excluding IMBHs within 8 
kpc of the VLII halo center. The resulting distribution of 



peak brightness for the IMBH stellar clusters would have 
to be further scaled down by a factor of 5, to account for 
the limited sky coverage of the SDSS. After accounting 
for these corrections, only a handful of clusters should 
be detectable with current surveys. This is consistent 
with a recent search for stellar clusters surrounding re- 
coiled IMBHs in the SDSS source catalog that resulted 
in an upper limit of 100 ca ndidates in the Milky Way 
halo (jO'Learv k, Loebl 12012ft . We have also estimated 
the proper motion of IMBHs on the sky. Having tan- 
gential velocities ranging from to 250 kms -1 as seen 
from an observer at 8 kpc from the VLII halo center, and 
distances of up to 100 kpc, their proper motions fall in 
the range between 0.1 and 10 milliarcsec year -1 . This is 
similar to measured proper motions of distant halo stars 
and may be currently detec table with the Hubble Space 
Telescope (Sohn ct al. 2012), and future missions such as 
Gaia. 
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Fig. 7. — Angular size of stellar clusters around Galactic "naked" 
IMBHs in the "direct collapse" model (see text for details). The 
angular extent is calculated for an observer at f^GC = 8 kpc from 
the center of the VLII halo. Dashed line: angular resolution limit 
of current ground- based surveys. 
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